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Abstract 

We consider the canonical quantization of an ordinary fluid. The resulting long- 
distance effective field theory is derivatively coupled, and therefore strongly coupled in 
the UV. The system however exhibits a number of peculiarities, associated with the 
vortex degrees of freedom. On the one hand, these have formally a vanishing strong- 
coupling energy scale, thus suggesting that the effective theory's regime of validity is 
vanishingly narrow. On the other hand, we prove an analog of Coleman's theorem, 
whereby the semiclassical vacuum has no quantum counterpart, thus suggesting that 
the vortex premature strong-coupling phenomenon stems from a bad identification of 
the ground state and of the perturbative degrees of freedom. Finally, vortices break the 
usual connection between short distances and high energies, thus potentially impairing 
the unitarity of the effective theory. 

1 Introduction 

Empirically, all fluids we know of undergo a phase transition when we lower the temperature. 
Either they freeze, or they transform into more exotic systems, like super-fluids or Fermi 
liquids. Why are there no ordinary fluids at zero temperature? For weakly coupled systems 
we have a microscopic understanding of the low-temperature macroscopic behavior at finite 
density [1] (we also understand why such systems at high temperatures exhibit hydrodynamic 
behavior.) However for strongly coupled ones, such as the so-called non- Fermi liquids, we 
don't. Of course a strongly coupled system could look like anything at long-distances — 
there is no in-principle preference for the ordinary fluid dynamics. Nevertheless, classical 
hydrodynamics is so common in nature at high temperatures, and as we will see, it can be 
defined purely in terms of low-energy degrees of freedom and symmetries like an ordinary 
QFT, that it is natural to ask whether there exist strongly-coupled quantum systems that 
at zero temperature and finite density do behave like ordinary fluids. 

It is tempting to conjecture that we know of no such systems simply because the corre- 
sponding quantum effective field theory would be inconsistent. As usual, the advantage of 
this viewpoint is that, as long as we allow for the most generic local dynamics involving all 
long-distance degrees of freedom and compatible with the symmetries, we are allowed to be 
completely agnostic about the microphysics yielding such macroscopic dynamics. Further 
motivation to investigate the consistency of the ordinary-fluid effective theory comes from 
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the recent interest in the so-called holographic liquids at low temperatures. There, one deals 
with strongly coupled systems at finite density via a dual description in terms of classical 
gravity, which can be thought of as providing the microphysics behind these liquid-like states. 
Characterizing the long-distance dynamics of such systems is non-trivial however (see e.g. [2] 
and references therein). Our approach allows us to make progress on the (perhaps modest) 
question: can some of these low-temperature, finite-density systems behave like ordinary 
fiuids? Partial indication that the answer may be 'yes' comes from the results of [3], but [2] 
argues that these cannot be consistently interpreted as a sign of hydrodynamic behavior. 

Without committing to any models for the microphysics, we will argue that: (i) the 
effective theory of an ordinary fiuid is not consistent; (ii) the effective theory of an ordinary 
fluid might be consistent after all; (in) there is no guarantee that such effective theory be 
unitary at low energies. Clearly, the matter deserves further study. 

Before starting our quantitative analysis which substantiates these claims, we conclude 
this introductory section with a few qualiflcations. First, we will only consider fluids at zero 
temperature. Thus, our results cannot be readily exported to the flnite-temperature case. 
Indeed we know that at high temperatures ordinary fluids abound in the real world (hence 
the 'ordinary'), and they exhibit no funny quantum effects at long distances like those we 
are going to discuss. In these cases we expect quantum effects to be overwhelmed by thermal 
ones. It would be interesting to validate this expectation quantitatively — a task which we 
leave for future work. 

Second, we are going to neglect dissipative effects throughout our paper. That is, we 
will consider perfect fluids only. The reason is that dissipative effects in hydrodynamics, 
like for instance those parameterized by viscosity and heat conduction, are associated with 
higher derivative corrections to the perfect fluid dynamics — see e.g. [1]. Therefore in the 
far infrared, that is for processes taking place on sufficiently long length scales, they can be 
safely ignored. Moreover, one may expect that the actual coefficients weighing these higher 
derivative corrections approach zero when the temperature is taken to zero. For instance in 
[2] it is argued that a finite viscosity at zero temperature is incompatible with hydrodynamic 
behavior. So, it is conceivable that by working at low enough temperatures and at long 
enough distances, one can make dissipation doubly negligible. 

Third, a crucial role in our analysis will be played by vortices. Precisely the existence 
of 'light' vortices is what distinguishes an ordinary fiuid from a superfiuid at the classical 
level. As long as one concentrates on the compressional modes — the sound waves — both 
systems obey hydrodynamics [1], and this holds at the non-linear, relativistic level as well 
[5]. However in a superfiuid the velocity field is irrotational, which implies that any vortex- 
like configuration will be singular at the center of the vortex, along a line, with the curl of the 
velocity field behaving like a delta-function peaked on this line. This means that from the 
viewpoint of the long distance/low energy effective field theory, the vortices are really UV- 
degrees of freedom, with finite energy per unit length (which is in fact mildly IR- diver gent.) 
So, for instance, one cannot form vortices by scattering phonons of very low energies — there 
is a gap, and as long as one works below the gap, the vortex degrees of freedom can be 
ignored On the other hand, in ordinary fiuids, one can build vortex configurations that 

close relative of the superfiuid vortex is the roton excitation, which is also gapped, and which can 
then also be neglected in the far infrared. 
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are arbitrarily mild, that is, that involve arbitrarily low momenta only. As a consequence, 
there is no gap in the energy one can store in a vortex. Vortices in an ordinary fluid are 
low-energy degrees of freedom, and they belong in the low-energy/long distance effective field 
theory together with the sound waves. In fact, we will see that, in a sense to be made precise 
below, the vanishing of the vortex gap is stronger than that of the sound wave gap — vortices 
are 'more massless' than sound waves. This will be the origin of all the quantum-mechanical 
peculiarities we will discuss, which are therefore absent for a superfiuid. 

2 The classical theory 

To begin with, let us review how classical hydrodynamics can be cast into a field theoretical 
language. We will adopt the viewpoint and notation of [5], to which we refer the reader 
for details. In particular, we will parameterize the fluid's configuration space by giving at 
time t the comoving ('Lagrangian') coordinates (p^ of each fluid element, as a function of the 
physical ('Eulerian') position x occupied by that fluid element: 

(f)' = (f)'{x,t) , 1 = 1,2,3. (1) 

Of course this description is completely equivalent to the inverse one, whereby one gives 
X as a function of the comoving coordinates and of time, and which we will also use at 
some point. However we find our starting point more convenient to construct the theory, 
because it identifies the fluid's macroscopic degrees of freedom with three scalar functions 
of spacetime coordinates, i.e. with three scalar fields: this way keeping track of Poincare 
invariance is straigthforward, and so is coupling the fluid to other systems, like gravity for 
instance. Moreover, as we will see in a moment, hydrodynamics follows straightforwardly via 
standard effective field theory (EFT) logic once we identify the correct internal symmetries. 

Before proceeding, it is worth stressing that we will be dealing with a fully relativistic 
theory, even though most laboratory fluids are highly non-relativistic. For these one could 
impose Galilean invariance rather than Poincare invariance, but this would not simplify the 
analysis we are going to carry out, conceptually or algebraically. We thus see no reason why 
not to keep track of relativistic effects and just neglect them when appropriate. 

We now come to the symmetries. The spacetime ones are of course the Poincare group, 
under which our 0^'s transform as scalars. As for the internal ones, we have a huge re- 
dundancy in choosing the fluid's comoving coordinates. This is not a symmetry — it is the 
standard arbitrariness one has in parameterizing a Lagrangian system's configuration space. 
To make any progress, we should make an explicit choice. A particularly convenient one is 
the following: At some given reference pressure we demand that for the homogeneous and 
static fluid configuration — the fluid's 'ground state' — the comoving coordinates be aligned 
with the physical ones: 

(f)' = x' . (2) 

It is then clear that homogeneity and isotropy for the physical properties of such a state 
cannot emerge unless the dynamics are invariant under internal translations and rotations: 

0^ -> + (3) 
^ O^0^ (4) 
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where O is an S'0(3) matrix. So far, we have not specified what distinguishes a fiuid from 
an isotropic sohd (a 'jelly'). It is an additional symmetry — the invariance under volume- 
preserving diffeomorphisms — 

0^ ^ , det = 1 . (5) 

This should not be confused with a trivial relabeling of the comoving coordinates, which we 
already got rid of. Rather, it corresponds to an invariance of the dynamics under physically 
moving fiuid elements around without compressing or dilating the fluid anywhere. If we 
were to do so in a solid, we would feel transverse stresses trying to pull all volume elements 
back to their rest position. In a fluid, on the other hand, we only feel reaction forces against 
compression or dilation. 

Invariance under shifts (eq. ([3|) forces each fleld (j)^ to be acted upon by at least one- 
derivative. At low momenta/low frequencies, the most relevant terms are those with the 
fewest derivatives. Therefore, the lowest order low-energy Lagrangian will involve exactly 
one derivative acting on each 0^. Poincare invariance then forces the Lagrangian to depend 
on the matrix 

= d^<P^ d^'cj)-' (6) 

only. Internal rotations (eq. Q) impose that we focus on 5*0(3) invariant functions of B^'\ 
and the volume preserving diffs (eq. ([s])) select the determinant among these. We thus have 
that the most generic low-energy Lagrangian compatible with all the symmetries is [5] 

S= [d^x F{B) , 5 = det 5" , (7) 



where F is a generic function. 

It is straightforward to check that the action ([T]) describes the dynamics of a perfect fluid. 
The stress energy tensor is 

T^, = -2F\B)B Byj d^<j)'d,<f)' + r]^.F{B) (8) 

(we are using the 'mostly plus' signature for the metric), which matches the standard form 
TfMu = (p + p)ui_iUi, + prj^iy upon the identiflcations [5] 

p = -F{B) , p = F{B) - 2F\B)B , = -^e^'^f'^ejjK d^cp'd^cp'd,^'' . (9) 

6v B 



In particular, we see that both p and p depend just on the degree of compression B, or 
equivalently, p depends on p only — our fluid is 'barotropic'. Different choices for F{B) thus 
correspond to different equations of state p{p), and once the equation of state is given, 
F{B) is uniquely determined. Notice that p, p, and u'^ are all invariant under our internal 
symmetries, eqs. ([3]-[5]), and so is T^^,. In fact, u^^ is invariant under generic internal diffs, 
with no volume-preserving restriction. What matters for characterizing the fluid flow is just 
that comoving coordinates do not change along it. Such a requirement is clearly preserved 
by generic diffeomorphisms of the comoving coordinates. 
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Since we have the correct stress-energy tensor for a fluid, we also have the correct hy- 
drodynamical equations, which follow from stress-energy conservation. The classical ground 
state of the fluid — the equilibrium configuration at a given pressure or density — is given by 
eq. ([2]). This spontaneously breaks all of our spacetime and internal symmetries, except for 
the diagonal combinations of internal shifts and spacial translations, and of internal rota- 
tions and spacial ones. As a result, there are gap less Goldstone bosons — the phonons — of 
which only the longitudinal one propagates. Indeed, we can study the propagation of small 
perturbations of the ground state by splitting (p^ = + and expanding the action at 
second order in the vr's. We get [5] 

S, = jd'x{-F'{l))[\i^'-\cl{V-^f] (10) 

where we defined the coefficient as 

2 _ 2F"{B)B + F'{B) 

and we stopped differentianting between internal indices and spacial ones, since they trans- 
form in the same way under the unbroken combination of internal rotations and spacial ones. 
In other words, from now on we should think of vf as a spacial vector field. We see from 
the quadratic action for 7? that only its longitudinal component has a gradient energy. The 
corresponding free solutions are plane waves propagating with speed Cs — the speed of sound. 
From the expressions for p and p as a function of B, eq. one realizes that = dp/dp\^^^, 
thus making contact with the usual expression for the sound speed in a perfect fiuid. For 
a non-relativistc fiuid <C 1, whereas for an ultra-relativistic one cj. ~ 1/3. We will not 
commit to either case, but instead leave cj. as a generic parameter. 

The transverse excitations do not have a gradient energy and as a consequence obey a 
free particle-like equation of motion, whose general solution is linear in time: 

7fy = V X (a(f) + b{x) ■ t) , (12) 

— * 

where a and h are arbitrary vector functions. This is the linearized limit of a vortex in 
constant rotation. For this reason we will refer to the transverse excitations as 'vortices'. 
Their lack of gradient energy is, of course, a direct consequence of the volume-preserving 
internal diff invariance, eq. ([s]), and is at the origin of all the peculiarities we are going to un- 
veil. Notice, however, that our diff-invariance is not a local symmetry, and as a consequence 
the configurations spanned by it — the vortices — are not gauge-modes, but real dynamical 
degrees of freedom. For instance, they have non-vanishing conjugate momenta; they just do 
not feature wave solutions. 

A more complete analysis of this classical field theory is carried out in ^ , where a num- 
ber of non-trivial results are derived — most notably a relativistic generalization of Kelvin's 
theorem, and the equivalence between the zero-vorticity sector of our fiuid and a superfiuid 
(i.e. a derivatively coupled scalar with a time-dependent background.) Here, instead, we will 
consider the quantum theory, and try to make sense of it. 



(11) 
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3 The naive effective theory 



The structure of the quadratic Lagrangian ( |10| aheady signals that, upon canonical quan- 
tization, we might be facing a strong-coupling problem for the vortices. The reason is the 
following: Consider first as a toy model a quantum-mechanical oscillator with some anhar- 
monic corrections to the potential. In perturbation theory, one first solves the harmonic 
problem, thus getting the standard oscillator spectrum, and then treats the anharmonicities 
as small corrections. The approximation is justified for those states whose wavefunctions 
are localized in a region where the potential is dominated by its quadratic approximation. 
So, for perturbation theory to be applicable in this case, one needs at least the ground state 
to have a localized enough wave-function (highly excited states will always be outside the 
regime of validity of perturbation theory.) Of course, what localizes the ground state is the 
curvature of the harmonic potential — the oscillator's frequency. For the system to be 'weakly 
coupled', one thus needs a steep enough quadratic potential. If we now move on to field the- 
ory, the role of the quadratic potential is usually played — in the absence of mass terms — by 
the gradient energy. For given spatial momentum k, the gradient energy gives a potential 
oc /c^lv^fcP- The vacuum wavefunction is thus localized about ipk = 0, and cannot probe large 
field values where interactions may become important. In the absence of a gradient energy, 
on the other hand, each mode's vacuum wavefunction is totally delocalized in the quadratic 
approximation, and its dynamics are completely determined by the interactions. We thus 
reach the conclusion that a (massless) field theory without gradient energies is prone to 
strong coupling, at all scales. 

There is a number of caveats in applying the above logic to our case. The first is that 
the absence of gradient energy may be an accidental feature of the lowest order in the 
derivative expansion. This is the case, for instance, for the ghost condensate [6], where 
gradient energy starts at the four-derivative order, -Egrad oc (V^vr)^. In the absence of 
quadratic terms with fewer spatial derivatives, such a term cannot be relegated to the class 
of higher-dimension operators, because it is marginal by definition — together with the kinetic 
energy i^kin oc tt"^ it determines how things behave under rescalings. In this case then, there 
is a well defined perturbative expansion. But this way out is not available to our vortices: 
the absence of gradient energy for them is enforced by a symmetry, which also forbids higher 
spatial-derivative quadratic terms. In the absence of time-dependence, exciting vortices 
costs nothing: we can deform the ground state (p^ = in the 'transverse' direction via 
eq. § and pay no energy price, and this extends to non-linear order as well. The second 
caveat, more relevant for us, is that the above quantum oscillator toy model assumes that 
the anharmonic interactions are of the potential form — only in this case delocalization of 
the wavefunction necessarily leads to strong coupling, because having access to large values 
of q entails having access to large interactions. But in our case, by construction, we only 
have derivative interactions, and moreover the very same symmetry that forbids the vortex 
gradient energy is also going to forbid many interactions involving vortices. In particular, 
as we will see more concretely in the following, all vortex interactions that do not involve 
at least two time derivatives are forbidden. Therefore the connection between wavefunction 
delocalization and strong-coupling is less obvious in our case. 

To settle the question, we should probe the theory by computing some physical quantity 
and check whether the perturbative expansion holds. The ideal candidates are usually 5*- 
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matrix elements, but here we face a complication. The longitudinal phonon has standard 
wave solutions, which upon canonical quantization, get mapped onto standard free-particle 
states. The transverse phonons, in contrast, do not behave as waves, and as a consequence 
there are no quantum asymptotic states associated with them. The classical field ttt behaves 
like a collection of infinitely many free particles rather than infinitely many oscillators. Upon 
quantization, its Hilbert space is not made up of standard Fock states. Without asymptotic 
states there is no S-matrix. 

A possible alternative, is to compute instead local n-point functions in real space, and 
to check whether perturbation theory holds for them. They may be as physical as the S- 
matrix: they characterize the physical interaction among local sources that couple to our 
fluid. We do not need asymptotic states to set up such a question. For instance, we can 
define the theory and the associated correlation functions via the path-integral formulation. 
Another possibility, which we will choose, is to give the theory asymptotic states for the 
vortex degrees of freedom by deforming it in the IR. We can add to the classical action a 
term that is compatible with all the symmetries except for the volume-preserving diffs, 

A^ = F'(1) jd^x\clB" , cl<^cl (13) 

and whose only effect, once expanded about the ground state, is to introduce a small gradient 
energy for the transverse Goldstones 

Jd^x{-F\l))[l^^-lcl{V-nLY-l4{V^7r'T'^,4)] . (15) 

We thus have wave solutions, propagating with speed ct, for the 'vortices' in the deformed 
theory — we promoted the vortices to real transverse phonons. Essentially, we are deforming 
the fluid into a solid/jelly that is stiff under compressional stresses but very soft under 
transverse ones. The original theory is recovered in the ct — >■ limit, with a qualification of 
course. With this ct deformation we are perturbing drastically the far infrared of the theory. 
We are adding asymptotic states, and we are going from not having an S'-matrix to having 
one. So from this viewpoint the fluid limit is obviously discontinuous. However, we expect 
more local quantities like n-point functions to be continuous in this limit. The situation 
should be similar to having a fairly narrow unstable particle: strictly speaking it is not an 
asymptotic state, yet for processes happening at time- and distance-scales much shorter than 
the particle's lifetime, we can treat it as an asymptotic state and associate an S'-matrix to it. 
For scattering processes faster than ~ 10 minutes, neutrons behave like asymptotic states. 

So, concretely, here is our program. We will consider scattering and decay processes in 
the CT-deformed theory. In particular, for simplicity we will stick to processes that involve 
at most four external legs. Thus to carry out calculations at tree-level, we need to expand 



^More precisely, the expansion of B^^ is 

= +2V • 7?+ (VjTT^' VjTT^) . (14) 

The hnear term is a total derivative, and can thus be neglected. The other terms, on top of giving the 



transverse phonons a gradient energy, correct the kinetic and gradient energies already present in ( 10 1 
However in the limit c^, <C < 1 these corrections are also negligible. 
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the action up to quartic order in the vr field; this is done in the Appendix, and the result 
is reported below. By construction, interactions involve one derivative per field. For finite 
ct, the theory is a standard derivatively coupled theory, and thus strongly-coupled in the 
UV. The strong coupling scale will depend on all parameters of the theory; however, we are 
interested in the c^-dependence, since eventually we will be taking ct to zero while keeping 
everything else fixed. If the strong coupling energy scale slides to zero in this limit, or 
equivalently, if cross sections and decay rates at fixed momentum or energy blow up in this 
limit, the theory is strongly-coupled at all scales, and thus inconsistent. Notice that we are 
trying to ascertain the consistency of the theory by computing something — the S'-matrix — 
that loses its meaning in the limit we are interested in. Still at finite ct we expect that 
the strong-coupling scale for the S'-matrix be related to a similar strong-coupling scale for 
n-point functions — that is, that for the latter perturbation theory break down at a distance- 
scale given by the strong-coupling scale inferred from the S'-matrix. Thus, even though the 
S- matrix does not exist in the ct — )■ limit, the formal fact that it appears to be strongly 
coupled at all scales is probably signaling that real-space n-point functions cannot be reliably 
computed in the fiuid theory, at any distance scales. 

It is convenient to rewrite the original Lagrangian ([T]) as 



+ H1 



£ = -wo/(V5) (16) 

where Wq = — 2F'(1) = {p+p)b=i is the ground-state's enthalpy density, and / is normalized 
accordingly, so that /'(I) = 1. With this new notation the speed of sound (11) is simply 
= /"(1)//'(1) = /"(I). Note that the derivatives here are with respect to V B. Also, we 
will use Svr to denote the matrix with entries {d'K)ij = diiij, and the brackets [ • • • ] to denote 
the trace of the matrix within. Then, up to fourth order the action is (see the Appendix) 

lldrr] [dn'] - | (3c^ + f,) [dnf + ^(1 + c^) [d7i]^' - if ■ dn ■ if 
. clldTi] det dn - Ic'Adn'Y + l{cl + fs) [dn'] [dnf - ^ (3c^ + 6/3 + U) [dn]' 
+ TT • dir^ • TT — (1 + c^)[(97r] if ■ On ■ n + ^\dn^ ■ 7f|^ 

+ + 3c^ + h) [dn]' - (1 + cl) [d7r'])7f' + |(1 - c^) if'] . (17) 

The first line is the free part of the Lagrangian, including the c^-deformation. The second 
line collects the trilinear interactions, whereas the third and fourth lines collect the quartic 
ones, /s and stand for /"'(I) and /""(I), respectively. Finally, notice that via the suffix T 
we indicate the transpose of a matrix, rather than the transverse part of tt as we did above. 

At this order we have four free parameters: c^, c^, /a, and f^. The dimensionful overall 
factor of Wq just gives us some reference units — we could use units in which it is one. For 
c^, we know that we want cj.. As to c^, it will be much smaller than one for a non- 

relativistic fiuid, and of order one (1/3) for an ultra-relativistic one. In the former case we 
expect /s and to be naturally of order c^, or smaller. If they were larger, B = 1 would be 
a special point for the shape of f{\/B), since by going to, say, 5 = 2, the second derivative 
of /, which controls c^, would undergo a relative change of more than order one. Likewise, 
in the ultra-relativistic case, for the same reason we probably want /a and of order one. 
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or smaller. That is, if we assume that 5 = 1 is a fairly generic point for /, and have 
to be at most of order c^. On the other hand it may be possible to have a fluid with some 
feature in the equation of state, where /" is small but higher derivatives of / are large. In 
the following we will make no assumptions about these couplings, since even carefully chosen 
values for them do not lead to drastic simplifications for our computations. Similarly, the 
non-relativistic case ^ 1 is only slightly simpler that the fully relativistic one, and we 
thus see no reason why not to investigate the latter. 



3.1 Sound-wave strong-coupling scale 

As a warmup, we estimate the strong-coupling scale for longitudinal phonon scatterings 
by ignoring the vortices, both as external states as well as internal lines. We assume for 
simplicity that /a and are of their 'natural' size, c^. We also assume that there are no 
cancellations among the various interactions. So, schematically the structure of the action 
is 

5,ound ~ J d^X dt Wo [{n^ - Cl d^TT^) + + 8^71^] , (19) 

where d stands for a typical spatial derivative, and for the interactions we used tt ~ Cs dn, 
valid for not terribly off-shell phonons. Now, to estimate the size of the amplitude at a given 
energy or momentum, we can proceed as follows. First, we redefine the time variable 

t t/cs , (20) 

to get a relativistic-looking kinetic term: 

S-sound ~ j d^x dt/cs Wo [cli-k^ " dV) + c] dV + SV^] (21) 

~ woCs j d^x [{-k^ - dV) + dV + . (22) 

Then, we notice that Cg has factored out of the action, and that combined with wo it gives 
the only energy/momentum scale in the action: = WqCs- The rest has standard rela- 
tivistic scaling (without Lorentz-invariant contractions though), which means that we can 
apply standard relativistic amplitude estimates. All interactions inside the integral have 
unit coefficient; the typical 2 — )■ 2 dimensionless amplitude is thus k'^, combined with the 
appropriate powers of M to match dimensionality 

interaction strength ~ ~ , (23) 

WqCs wocl 



■^This is what happens for instance for the constant (i.e., _B-independent) Lagrangian 



C = ^woiVBy^^ , (18) 

which corresponds to the simple equation of state p = cj.p. Notice however that such a simple case, besides 
being extensively considered by cosmologists, is not preferred in any sense over more generic equations of 
state — with the exception of the ultra-relativistic case, where the linear equation of state p = 1/3 p follows 
from scale-invariance. 
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where E is the typical energy in units of the original time variable (notice that spatial coor- 
dinates are untouched, so there is no such ambiguity for k.) The strong-coupling momentum 
and energy are thus 

/c* = [wqCsY^'^ , E^ = Csk^ , (24) 

respectively. 

The above estimate yields the correct strong-coupling scale for longitudinal phonons. We 
could do the same for the vortex sector, at finite ct- However as we will see, for the vortex 



interactions there are cancellations that are not manifest in eq. (17) and that would impair 
this simple estimate. 



3.2 Hunting for all factors of c 

When we start computing amplitudes and physical quantities like cross sections and decay 
rates, we have to be careful about extra factors of Cs and of ct besides those appearing 
explicitly in the various interaction terms. For instance, we just saw that the longitudinal 
phonon interaction strength is of order k'^/{wQCs), whereas the Lagrangian interaction terms 
are proportional to c^. In hindsight, this result just follows from dimensional analysis, once 
we keep separate units for space and time (we can still set h = 1 though.) The quartic 
interactions involve four powers of momentum. The interaction strength is thus fc^ divided 
by whatever combination of wq and has the same units as k^. Of course, we have an 
ambiguity as to the units of Wq — is it a mass- or energy-density? We can easily resolve this 
ambiguity by looking at the kinetic term. By construction our vr^ have units of length; the 
action is dimensionless (for h = 1); wq is thus a mass density, and wqcI an energy density: 
[poc^] = Ek^ = Csk^. The dimensionless combination therefore is Csk^/{pocj.) = k^/{poCs), as 
expected. 

So, a possible strategy to get all the factors of Cg right in amplitudes, cross sections, and 
rates, is to use the standard relativistic formulae, and then insert suitable powers of to 
match dimensions. Essentially this is equivalent to redefining the time variable as we did 
above, to end up with a relativistic kinetic term with Cg = I. However this strategy is going 
to fail once we include vortices/transverse phonons in our processes: with two different prop- 
agation speeds and ct for longitudinal and transverse polarizations, dimensional analysis 
does not suffice. Equivalently, by redefining the time variable we can cast only one of the 
two kinetic terms in relativistic form. 

In the Appendix we will therefore briefly review the standard relativistic formulae and 
derive the modifications needed to apply them to our case with c 7^ 1. We adopt this 
somewhat cumbersome action plan, rather than going through some standard condensed 
matter textbook and trying to dig up the relevant non-relativistc formulae, for no other 
reason than we are more familiar with the relativistic Feynman rules and related formulae — 
and we assume that the reader is also. The bottom line is pleasantly surprising: We can 
use the standard relativistic Feynman rules and formulae for infinitesimal cross-sections and 
rates, with no modifications, even when we start considering different fields with different 
speeds. By 'standard relativistic rules and formulae' we mean those associated with the 
so-called relativistic normalization of single-particle states, as derived for instance in Peskin- 
Shroeder [7]. 
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As to the overall factor of Wq in (17), it is straightforwardly kept track of. Either by 
inserting for each Feynman diagram a wq for each vertex, a. 1/wo for each internal line, and 
a 1/ y^wo for each external line, or most simply by setting it to one and retrieving it at the 
end of the computation via dimensional analysis. 



4 Simple processes — the vortex strong coupling 

We compute, to tree level, a number of simple processes in order of increasing number of 
vortices on the external legs. For scattering processes, for simplicity we will only consider 
initial states with zero total momentum. Given that Lorentz boosts are spontaneously broken 
and that we have a preferred reference frame, this is a non-trivial choice — we are setting 
some kinematic invariants to zero. With an abuse of language, we will refer to this choice 
as "working in the center of mass (CM) frame." For the decay of a single finite energy 
excitation, on the other hand, this choice is not an option, of course. We will use the 
formulae for amplitudes, cross sections and rates found in the Appendix. But, as commented 
on before, the good news is that these formula look just like the usual relativistic ones that 
we are used to dealing with. So, except for the additional factors of Cg and ct coming from 
the external states' dispersion relations and from the internal lines' propagators, everything 
goes through just as usual: each external line carries a polarization-vector factor (times 
I/^/wq), each incoming or outgoing time-derivative contributes a ^iu, each incoming or 



outgoing spacial-gradient contributes a ±ik, and so on. From eq. (17), we immediately get 
the Feynman propagator: 



{Tn\x)n'iy)) - ■ , ^ . + - ■ , , (25) 

Wo uj^ — cjp^ + le Wo — c^p^ + ze 



where Pl and are the longitudinal and transverse projectors, respectively. 

We will not content ourselves with amplitudes. Rather, we will compute physical, mea- 
surable quantities like cross sections and decay rates. The reason is that amplitudes depend 
crucially on the normalization chosen for the single-particle states. For instance going from 
the so-called relativistic normalization to the non-relativistic one, would move some factors 
of Cg and ct from the amplitudes to the phase-space elements, in such a way as to keep cross- 
sections and rates unaffected. Ascertaining the strong-coupling of the theory in the — 
limit at the level of amplitudes requires a derivation of partial waves, a la Jacob- Wick, being 
careful about the factors of Cg and ct- Although we have also derived our results using that 
method, we found it simpler to present them by focussing on cross-sections and decay rates. 

A final remark about external vortices. When we take the ct — )■ limit we have to decide 
whether we are going to keep their momenta or their energies fixed. The first choice is the 
more conservative, since it corresponds to taking their energies to zero, thus weakening any 
possible strong-coupling phenomenon we are going to encounter. It is also the only consistent 
one, since the alternative one would send the vortex momenta to infinity, outside the regime 
of validity of any effective theory. In the following we parameterize everything in terms of 
momenta rather than energies, so that taking the ct ^ limit is straightforward. Notice 
also that only if we keep the vortex momenta fixed is our deformed theory with small ct close 
to the fluid one with ct = 0: in the Lagrangian weighs the gradient energy, so that by 
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sending ct to zero while keeping the momenta fixed one is in fact sending the magnitude of 
that Lagrangian term to zero. Related to this, it is somewhat tricky to deal with processes 
that include longitudinal phonons in the initial state but no longitudinal phonon in the 
final state: the initial longitudinal phonons' finite energy should be divided among the final 
vortices, thus making their momenta diverge for Ct 0. In other words, one cannot send 
Ct to zero and keep all momenta fixed. We will see an example of this below, in the decay of 
a longitudinal phonon into two transverse vortices. We postpone a discussion of the related 
subtleties until then. 

We will use p 's to denote the momenta of the longitudinal modes, and k's and e's to denote 
the momenta and polarizations of the transverse modes. Our e's are real, thus corresponding 
to linear polarizations, and normalized to one (hence the 'hat'.) For longitudinal phonons 
the polarization vector is p, of course. For all the processes we will just compute the leading 
contribution in the limit ct/cs ^ 1, for which we hope to learn something about the original 
fiuid (ct = 0). 



4.1 Longitudinal 2 — > 2 scattering 

This is the simplest of the scattering processes. To tree level, the only relevant diagrams are: 








We designate, here and for the rest of the paper, the solid lines as longitudinal excitations 
and the curly lines as the transverse excitations. Time fiows to the right. 

When done in the center of mass frame, the only kinematic variables are the momentum 
of the longitudinal phonons p and the scattering angle 9. To tree level, the total amplitude 
is given by 



iM 



LL~^LL 



■ A 2 

ip 

Wo 



[U/cl - 2fi/ct + 3cl + 2/3 + ci + 2(1 - 3c^) cos^ 



(26) 



Remarkably, the graphs with transverse propagators do not contribute to the amplitude, 
even individually. The infinitesimal cross section is 



da 



1 \M 



LL^LL\ 



cl 647r2(2p)2 



-dVL 



(27) 



where we made use of the phase space element computed in the Appendix (eq. ( |113 )). We 
can easily calculate the total cross section. The final particles are identical, so we over-count 
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when we integrate over all final phase space. To counteract this we simply include a 1/2 
symmetry factor. To all orders in Cs the total cross-section is 



2567rp2 \woc 



p 



+ — - + 



p 



P^ \WoCs 



(28) 



where a ^ [fjc] - 2fl/4 + Sc^ + 2/3 + 4) = 0(1) + ^(c^) + O(c^) (assuming h, h ~ c]) 
and/3 = 2(l-3c2). 

This result matches our dimensional estimate of the strong coupling scale for longitudinal 



phonons in eq. (24): the cross-section (28) is the geometric cross-sectional area for wave- 



packets of wavelength 1/p, times the square of the dimensionless interaction strength we 
estimated in sect. 3.1 We have strong coupling when a becomes of order l/p^ — in such 
a case the two wave-packets have an 0{1) probability of interacting — in agreement with 
sect. 13.11 



4.2 Longitudinal decay and vorticity production 

In addition to scattering cross sections, we can also calculate the decay rate of a longitudinal 
phonon into a longitudinal phonon and a transverse one. This is kinematically allowed thanks 
to the difference in propagation speeds between longitudinal and transverse excitations. To 
tree level, we simply have one diagram, the longitudinal-longitudinal-transverse vertex: 



Imposing the kinematical constraints coming from momentum and energy conservation, ex- 
panding in ct/csi and keeping only to first order in this parameter, we can write the amplitude 

as ^ 

iMl-,lt = Z^^I^^e ■ p) sin 6/2 (cos 6 - c]) (29) 



where is the angle between the outgoing longitudinal mode and the decaying one (with 6* = 
corresponding to forward decay.) If the vortex is polarized orthogonally to the scattering 
plane the decay amplitude vanishes because of parity conservation, whereas for parallel 
polarization it further simplifies to 

iMl^lt, = sin e (cos e - c^) (30) 



/Wo 

Notice the fairly featureless angular dependence, which combined with a similarly featureless 



phase space (eq. (116)) yields the surprising result that with an order-one relative proba- 
bility the longitudinal phonon will recoil backwards, i.e. with ^ > 7r/2, by emitting a high- 
momentum forward vortex. Moreover, the unusual kinematical constraints associated with 
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the ct ^ Cs hierarchy force the vortex to be emitted always in the forward half-space, and 
at an obtuse angle relative to the final longitudinal phonon. 



The physical quantity that we want to calculate is the decay rate F, given by (97). Using 



the infinitesimal phase space for a two-particle final state with non-zero total momentum 



and Ct ^ Cg, given by (116), summing over possible transverse polarizations and integrating 



over the solid angle we get: 

Ct 4 



L^LT 



Cs Wo SIStt 



(5-6c2 + 21c^) (31) 



which, as we can see, smoothly goes to zero as ct — )■ 0. 

In the same way we can study the decay L — )■ TT of a longitudinal phonon in two 
vortices. Notice that the decay kinematics in the ct/cs <^ 1 limit implies ki ^ —k2 and 
2\ki\ ~ \p\cs/ct ^ \p\- That is approximately two back-to-back vortices, carrying each half 
the energy and with momentum scaled up by a factor ~ Cs/ct with respect to the initial 
state phonon. We find the amplitude 

2 3 

iMl^tt = [ei ■ €2(1 + 2 cos Q'') ^ p ■ e,p ■ €2] (32) 



where represent the polarizations of the two vortices and cos 9 = p- ki. For the decay rate 
we have thus roughly 



p^cl 

L^TT ~ 5- 

WqC^ 



(33) 



corresponding to a 'quality factor' Vl_^tt/^ ~ P^Cg/woc|.. A phonon with momentum 

~ {woCsf/\cT/csf'' (34) 

has a width comparable to energy. We can thus identify strong interaction scale 

for longitudinal phonons. Notice that this scale is parametrically smaller than the naive 



estimate in eq. (24). This fact is largely a consequence of the peculiar kinematics of the 
decay L — )■ TT, where starting from an initial quantum of momentum p, the final state 
quanta have a much higher momentum scale k ~ pCs/cx- This higher scale is naturally 
associated with a stronger interaction strength. The fact that starting with soft quanta 
one can probe much shorter distances due to the large final state momentum also suggests 
more care with the use of the notion of effective field theory. We will elaborate briefly 
on this in sect. [6j Notice also that the corresponding vortex momentum strong scale is 
instead fc* ~ p*Cs/ct ~ {wqCtY^^{cs/ ctY^'^ . The computation of TT — )■ TT scattering in 
the next section will show that the vortex momentum cut off is actually ~ {wqCtY^^^ which 
is parametrically smaller. Then, within the resulting smaller range of validy of the effective 
field theory, the process L — )■ TT remains weakly coupled. 

The results just derived display one general property of amplitudes involving vortices: 
they are accompanied by at least one power of the vortex energy. |^ That property directly 



^In the case of the L LT decay this also corresponds to a significant suppression of the ampHtude, 
given the parametrically suppressed value of the vortex energy in that process. In the case of L — ^ TT the 
energy of L and T modes is instead comparable. 
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follows from invariance under volume preserving diffeormorphisms and can be made evident 
by chosing a suitable field parametrization. In a neightbourhood of = x the most general 
field configuration can indeed be written implicitly as 



(j){x,t) = g{y,t) y = x + 7rL{x,t) 



(35) 



where til = Vip{x,t) is a longitudinal perturbation while g{y,t) is a volume-preserving 
diffeormorphism generated by "exponentiating" a transverse vector field 7fjn{y,t) (Vy ■ ttt = 
0). That is 



lim 

Af-s>oo 



e + 



(f,t)=y + 7fr(y,t) + 0(4) (36) 



where e{y) = y is the identity function and o represent function composition. Our procedure 
to define a finite transformation g starting from the infinitesimal one y + ttt is just the 
standard exponential map of Lie groups. In the end, eq. (35) is a three dimensional (in 



field space) family of field configurations that for tt^ ^ — )■ reduces to the most general one 



X + ttl + ttt- It follows that in a neighbourhood of the identity, eq. (35) is a faithful (one to 



one) parametrization, acceptable to perform perturbation theory. Now a time independent 
7rT{y,t) = 7fT{y,0) is just a symmetry tranformation of the action in the limit Ct = 0. 
Therefore, apart from the 0{c^) kinetic perturbation, ttt enters the lagrangian with at least 
one time derivative, and therefore amplitudes have the corresponding suppression. Notice 
that in the computation of L — )■ LT, where we used the simple parametrization = x + vf , 
that result arose via a non trivial cancellation of different terms in the amplitude. 

As a further check of the above general property of amplitudes involving transverse modes, 
consider the process LL — )■ LT. The relevant Feynman diagrams are 








The amplitude vanishes if the vortex is polarized orthogonally to the scattering plane, 
whereas for parallel polarization it is 



iM 



LL^LTh 



Wo 



cot e [3 - h/cl - 2cl -h + 3c^ - 2(3 - c]) cos' 6] + 0{4) , (37) 



again vanishing, as expected, as ct — with p fixed. This result corresponds to a cross- 
section scaling as ct , 

.4 \ 2 



o'll^lt 



Ct I f p 



Cs \WoCs 
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(38) 



The results just derived deserve one additional comment. For a classical inviscid fluid, we 
know that if we start with zero vorticity everywhere in space, we cannot produce any. How- 
ever this is not because of some global charge conservation, but, more prosaically, because 
the source of vorticity is proportional to vorticity itself. For a non-relativistic fluid: 

n = -{v-v)n+ {n-v)v-n{v -v) . (39) 

and X, with mutual interactions 



This is similar to, say, a scalar field theory with two fields 
of the form the x e.o.m., the source term is proportional to x itself. 



6C 

Sx 



3 



(40) 



which means that x = is a perfectly good classical solution no matter what (j) does. On the 
other hand, we know that this fact does not survive quantum mechanically. Two quanta 
in the initial state can annihilate in a vertex to yield two x quanta in the final state. 
However this way one will never produce a single x quantum if this is not already present 
in the initial state. Coming back to our fluid, it is suggestive to interpret our results above 
in light of this analogy. We expect that the classical non-generation of vorticity will not 
survive at the quantum-mechanical level. Vortex excitations will be generically produced in 
scattering and decay processes, even if there are no vortices in the initial state. However, 
the production of a single vortex quantum should be prohibited, being immune from the 
aforementioned effect'. Indeed consistently with this expectation, at fixed phonon 



momentum p, we found the single vortex amplitude Eq. (29) vanishes for ct 0, while the 



two vortex amplitude Eq. ( 32 ) does not 



4.3 Longitudinal and transverse scattering 

Another interaction we can consider is the scattering of a longitudinal excitation and a 
transverse excitation. The tree level diagrams are given schematically by: 








Besides the usual scattering angle, there are also the additional degrees of freedom associated 
with the polarization of the transverse modes. The amplitude is 

i2csCTP'^ 



iM 



LT^LT 



Wo 



(41) 



^In fact, associated with the invariance under volume-preserving diffs there are infinitely many conserved 



local currents and global charges [F. Eq. (39) is a consequence of these infinitely many conservation laws, 



but it does not take the form of a simple conservation law for vorticity itself. 
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Note that M. oc ct, consistent with the general argument presented in the previous section. 
It should also be noted that there is explicitly no dependence on and /s. This result can 
be easily understood using the field parametrization discussed in the previous section. The 
terms proportional to and /s only depend on the longitudinal field ttl- Then obviously 
does not contribute to LT — )■ LT while /s can only contribute via the second diagram in the 
upper line in the figure, which is clearly of order c^. Again, in the standard parametrization 
all these results follow from a non-trivial set of cancellation, which also represent a check of 
our computations. 

The infinitesimal cross-section is given by 

""-IJIJ: (^) ( 4.u,o(c! + c,) )^{[<P-^''---l(^.-P')fe-P)}-^» (42) 

initial e final e ^ ^ V\ 1 < LJJ 

Here we are averaging over the incoming polarizations (hence the 1/2) and summing over 
the final ones. A good basis to do this in would be parallel to the scattering plane and 
perpendicular to the scattering plane. As we can see from the form of A^, any perpendicular 
component of the polarization does not contribute to the amplitude. Putting everything 
together, keeping all powers of Cs, we have the total cross section: 



V 



4 \ 2 



We now see the importance of M.lt^lt oc ct- This dependence on the "transverse speed of 
sound" is necessary to avoid a divergent physical quantity (here the scattering cross-section) 
as Ct — 0. Remarkably, the strong-coupling scale for this process is the same as for purely 
longitudinal scattering — cf. eq. (28). Note also the independence of the cross section on /a 
and /4. This cross section is generic for all fiuid types, regardless of the particular functional 
form of f{\fB). Its dependence on the particular fiuid model only comes through the speed 
of sound. 



4.4 Transverse 2 — > 2 scattering 



As we can see from formula ([94j) for the differential scattering cross section, the more trans- 
verse incoming and outgoing states the more a rate could possibly diverge as ct — 0. The 
results so far agree with that expectation. The only problematic quantity we encoutered is 
the rate for L — >■ TT, which has however a peculiar (singular) kinematics as ct — >■ 0. The 
processes with smooth kinematics were instead found to have a well behaved rate. In the 
T + L —7- T + L cross section we picked up a c^^ from the phase space of the outgoing 
excitation and we picked up a c^^ from one of the l/2ii^ normalization factors. It was thus 
critical that |A^lt->lt| oc in order that the cross section be well defined. 



From eq. (94) we can see that for transverse 2 to 2 scattering ^ oc Cj, |A^| . In order 



that our physical process be finite in the ct ^ limit we need (at least) oc c^. We 
find that this is not the case. In fact, simple power counting using the parametrization 
discussed in section 4.2 indicates the amplitude will be suppressed (at least) as c^. This 
expectation is confirmed working in the standard parametrization: by dramatic cancelations 
among Feynman diagrams the zeroth and first order terms vanish, and the leading non-trivial 
term is of order c^. The necessary tree-level Feynman diagrams are: 



17 




»OoooO« 






In the CM frame with || and _L denoting polarizations parallel and perpendicular to the 
scattering plane, respectively, at lowest order in ct the amplitude is 



iM 



TT^TT 



Wo 



X 



COS 26 for 
I ( cos 9 — cos 29) for 



(44) 



and zero for all other combinations of polarizations. In the J-HX case, 9 is the angle 
between the two ||-polarized phonons. Note that, as in the previous physical process, there 
is no dependence on and f^. So, once again, this process is generic for all fluid models 
regardless of the details of f{\^). 

After squaring the amplitude, we average over incoming polarizations and sum over the 
final ones. We have: 



7 E E 



Wq 



[2+ I COS 2^+ I COS 4^] , 



(45) 



initial e final e 

and so the total cross section, including a 1/2 symmetry factor, is 

i I L6 \ i / \^ 

U TT—^TT 



/13 



2567r Vis 



/ F \WqCtJ 



(46) 



As we can see, it blows up as we take cj^ — )■ indicating that our theory is breaking down. 
We emphasize once again the absence of the free parameters /s and in this result, which 
implies that the pathology just unveiled cannot be avoided by a judicious choice of their 
values. 



5 The infrared situation 

Our S'-matrix analysis indicates that the transverse degrees of freedom are strongly-coupled 
at arbitrarily low energies. However the strong-coupling phenomenon we unveiled is quite 
peculiar. We deformed the theory in the far IR by introducing a small deformation parameter 
Ct- This changes the asymptotic states of the theory, and we discovered that some of these 
get strongly coupled in the UV, at an energy scale that drops to zero when we recover 
the original theory — the limit — )■ 0. A vanishing ultraviolet strong-coupling energy scale 
suggests that our problem is probably more properly thought of as an infrared one — we may 
be approaching strong-coupling from the wrong side! That is, in the deformed theory at 
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finite ct we encounter strong coupling in moving to high energies, but since the deformed 
theory differs from the original one at low energies, it may be that the strong coupling scale 
is in fact a divide between the two theories — that there is no regime where the two theories 
look alike. If we stick to the original description, without ever introducing cr, we may realize 
that we have some form of strong coupling in the IR. 

The distinction we are putting forward may sound like a matter of definition, but it is not. 
A theory that becomes strongly coupled in the UV is simply not defined at energies of order 
of the strong-coupling scale and above — it needs infinitely many parameters for its definition. 
If this were the case for us, our theory would not be consistent, in any energy range. On 
the other hand, there are a number of ways in which perturbation theory can break down 
in the IR without impairing the consistency of a theory. There is for instance real QCD- 
like strong coupling, where perturbation theory does break down but the non-perturbative 
theory is perfectly well defined — it is just hard to solve! Or there are QED-like infrared 
divergences, which can be tamed by focusing on suitable infrared-safe observables, for which 
the perturbative expansion applies. Or there may be huge quantum IR fluctuations without 
necessarily implying large interactions, like for instance for would-be Goldstone bosons in 
1 + 1 dimensions [8j. This would signal a bad identification of the theory's vacuum state. 
And in general, to ascertain the consistency of an IR-problematic theory, at least in certain 
energy and momentum regimes, one can just put the system in a finite-size 'box' and consider 
its time evolution for a short time, in which case perturbation theory typically does not even 
break down. 

The first signal that our problem with the transverse degrees of freedom may be infrared 
in nature comes from considering quantum fluctuations about the semiclassical vacuum state 
with (0^) = . An order parameter that conveniently quantifies the amount of spontaneous 
symmetry breaking in a manifestly translationally invariant fashion is 

{d,<P') = 5; ~ 1 . (47) 

It is straightforward to estimate quantum fluctuations in this quantity. We decompose the 
fields as cj)^ = + vr^, and from the vr^ propagators, 

{Tn\x)^\y)) ^ ^ • , + - ■ , (48) 

where Pl and Pt are the longitudinal and transverse projectors, we get 

1 5 

^ ^ ' Wo to ^ ^ 

(doTT^dy) ~ — / (50) 

Wo 

{doTi^doTi^) ~ — p^w (51) 

Wo 

These are real-space correlators, and in the right-hand sides the dimensionless, order-one 
part of the Fourier transform, J dVt d\ogp dloguj e^^'"\ is understood. Also we are considering 
considerably off-shell (a;,p) pairs, by taking for instance the separation in real space to be 
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space-like with the respect to the sound speed (i.e., by working in Euchdean space.) The 



correlators (50, 51) behave in essentially the same way as for standard field theories in 



dimensions higher than 1 + 1: quantum fiuctuations in order parameters are damped at low 
momenta and low energies, and as a consequence in the IR there is spontaneous symmetry 
breaking. However the correlator (49) ruins this familiar picture: for fixed momentum, it is 
IR-divergent. Equivalently, in the p-u plane there is a sector extending all the way to p = 0, 
u = where quantum fiuctuations in our order parameter (47) are huge. This reminds us of 
Coleman's theorem p] in 1 + 1 dimensions, and suggests that our strong-coupling problems 
may stem from an incorrect definition of the vacuum. In other words, we have been assuming 
all along that doing perturbation theory about the semiclassical vacuum 0^ = is sensible, 
but eq. (49) is somehow telling us that quantum fiuctuations want to dismantle this state. 
As an aside, notice that we have all reasons to believe that Lorentz boosts are spontaneously 
broken: since the correlators ( [50| [5l| ) are damped in the IR, it seems that only the breaking 
of the spatial symmetries (translations, rotations, volume-preserving diffs) is affected by our 
phenomenon. More concretely, we can compute quantum fiuctuations in an order parameter 
that is invariant under all symmetries but Lorentz boosts, like for instance the fiuid velocity 
u^. Classically in the ground state we have = (1,0). At first order in the vr field the 
fiuctuation reads 6u'^ = (0, — vr). From eq. (51) we thus get the velocity- velocity correlator, 

1 



~ 5f —p^u , (52) 



Wo 

which is damped in the IR, thus signaling that quantum fiuctuations do not want to restore 
Lorentz invariance. 

Several natural questions arise: How do we check the above statements more concretely? 
Can we identify the correct vacuum state? Does it support a well defined perturbative 
theory? And perhaps the most physically relevant: is there a semiclassical limit where we 
recover classical hydrodynamics? To start addressing these questions, we step back from 
our fiuid case and consider a much simpler system with somewhat similar features: the free 
particle in quantum mechanics. 



5.1 The free quantum particle 

Consider a free particle living on a line. Classically, x = Xi = const is a perfectly good state. 
However, quantum mechanically we know that if we start from a state localized around 
Xi, time-evolution will make the wave-function spread out, and at very late times the state 
will be totally delocalized. Related to this, the ground state of the theory wants to have a 
constant wave-function throughout the whole line. (This is non-normalizable for an infinite 
line, but for simplicity we can replace the line by a very large circle.) The spontaneous 
breaking of translations that we see classically, x = const, quantum mechanically is gone. 

Given the system's simplicity, there are many equivalent ways to describe this phe- 
nomenon quantitatively. One that will prove readily exportable to the fiuid case, is the 
path integral one. Consider, in the path-integral representation, the amplitude for evolving 
from Xi at t = to X f at t = T: 

{Xf, T \x„ 0) = / Vx e^^[-] oc exp ^^"^^^ . (53) 

Jx{0)=Xi ^-^ 
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We are using h = 1 units. Also for simplicity we are setting the particle's mass to one. The 
exponent is of course the classical action for the solution interpolating between Xi and Xf 
in time T. The prefactor that we omitted is the determinant of the kinetic operator about 
the classical solution. Since we have a free theory such an object does not depend on the 
solution under consideration, or on the overall distance {xf — Xi), and will play no role in 
our discussion — it is just an overall normalization factor. 

We see that, for a given T, at small distances \xf — Xi\ < \/T the amplitude is essentially 
constant, whereas it starts oscillating rapidly when we move to much larger distances, \xf — 
Xi\ ^ y/r. We are tempted to conclude from this that the width of the wave- function 
grows like a/T, but this is premature because the amplitude we got is a pure phase at all 
distances, which means that the probability is uniform in Xf. In other words, as soon as T 
is bigger than zero, the wave function is completely delocalized. This is clearly not what 
we expect physically, and is an artifact of having chosen an initial 5-like wave-function: a 
S{x) has infinite momentum spread, which means that in zero time the particle is going to 
be everywhere with non- vanishing probability. If we instead choose a more reasonable state, 
say a gaussian centered at = with width a, the final wave function we get is: 



,.2 



exp-2^ , T < 

'ip{xf,T)= I dxiip{xi,0) {xf,T\xi,0) (X { (54) 

exp — ^ , T a'^ 

where we ignored overall normalization constants, as well as factors that do depend on 
Xf but are pure phases. We thus see that at early times the wave- function width is — not 
surprisingly — dominated by the initial width a, whereas at late times it grows linearly with 
T. 

We can turn the problem around and ask: suppose we want to have a somewhat localized 
state that looks stationary over a time period of order T. What is the minimum wave-function 
width we should allow for? Answer: \/T. As a consequence, in the long-time limit there 
cannot be spontaneous symmetry breaking. With hindsight, we can go back to our original 



amplitude eq. (53) and learn how to read it properly. Despite being purely imaginary, the 
exponent really tells us the minimum uncertainty we should allow in the initial position in 
order for our state to look approximately stationary over a time of order T. 

Before moving on to the field theory case, it is worth pointing out that if we have 
free particles we can run the above computation independently for each of them, since all 
amplitudes factorize. This means that, roughly speaking, each particle can afford an 0{1) 
action, so that the whole system can explore trajectories with an 0{N) action. 



5.2 Coleman's theorem from the path integral 

We can now apply the same logic to a massless free field theory in d + 1 dimensions, 

S[<p] = jd^+^x\{d<P)\ (55) 

and recover Coleman's theorem in the path-integral language. We have a shift symmetry 
— )■ + c, which classically is spontaneously broken by any Poincare invariant configuration 



21 



= const. Without loss of generality we can set the constant to zero, and ask whether 
quantum fluctuations tend to disrupt this configuration. Coleman teaches us that this will 
be the case for d = 1. 

We want to compute the path integral 

(/(f),T|0,0) = / P0e*^[<^l oc e^^l-^^'l , (56) 

for a generic final configuration f{x). Given what we learned for the free particle on a line, 
the result will tell us which are the field configurations that inevitably get populated after 
a time T if we start at t = with a wave-functional for (j){x) centered around = 0. To 
compute ^[^ci], we have to solve the classical equations of motion with the given boundary 
conditions. In (spacial) Fourier space the immediate solution is: 

'/>ci(fc,^) = /(^)-— (57) 

where f{k) is the final configuration's Fourier transform. The action is 

S[(P,i]= jd'^k\f{k)\^ -k cot kT . (58) 

We can for simplicity consider the two regions /c ^ 1/T and k ^ 1/T separately, and ignore 
what happens at intermediate momenta. As to the former region, we can set cot kT — ?■ 1/ kT, 
so that we have 

S{k<^l/T)^^j d''k\f{k)\\ (59) 

The high momentum region, however, contains all the poles of the cotangent. These poles 
are just enforcing the right periodicity for the corresponding momenta. That is, for a given 
T all modes with k = rni/T should go back to zero after a time T, as enforced by the perfect 
harmonicity of our system. As a consequence, if the final f{k) does not vanish for these 
special momenta, the corresponding action is infinite. To have a finite action f{k) has to 
have infinitely many zeroes, at the right locations. Apart from this peculiarity, when we 
integrate over momenta much larger than 1/T the cotangent behaves like a number of order 
one, and to get an estimate for the action we can just ignore it. We thus have 

S(A;>1/T)~/ d'^k\f{k)\^ -k . (60) 

Jl/T 

We can now consider a specific final configuration f{x). We choose it to be localized 
in a region of size T, with Fourier momenta of order 1/L and with typical magnitude /. 



As to the final remark of sect. 5.1, here we are considering just a few independent degrees 



of freedom — the Fourier modes with k ~ 1/L m volume of size L — so that they get 



spontaneously excited to the desired level / only if the corresponding action in (56) is of 
order one, or smaller. At early and late times we get, respectively 

S{T « L) ~ , S{T » L) ~ L^-^p (61) 
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Our path- integral formula ( 56 ) at late times thus yields 



(/(f),T|0,0) -expiL'^-V^ T>L (62) 

This is Coleman's theorem in our language. For d > 1, the amplitude for finding the system 
significantly away from = (large /) in bigger and bigger regions becomes smaller and 
smaller. Put another way, the typical spread of the field decreases with distance scale. This is 
the standard behavior of a field theory in high dimensionality. On the other hand, for d = 1 
the L-dependence drops out, and we are left with a finite amplitude of finding an order-one 
/ on all scales. In fact we know that a more careful estimate shows that the typical spread 
of the field grows logarithmically with distance scale — as usual one cannot get the logs right 
by performing Fourier transforms by dimensional analysis. Our analysis for the fluid will be 
insensitive to such minutiae. 

5.3 Back to the fluid 

For our fluid, it is easy to convince oneself that the situation is more similar to the free particle 
case than to a 1 + 1 standard field theory. The reason is simply that all the peculiarities 
associated with the transverse excitations in a perfect fluid stem from their not having a 
gradient energy, that is, from their behaving like infinitely many free particles. Of course, this 
is a linearized statement. At the non-linear level, for large vorticous excitations, interactions 
become important and the free-particle approximation breaks down. But, in fact, it is easy to 
see that at low fluid velocities we can really describe the dynamics of vortices as a system of 
free particles constrained to move on a quite non-trivial infinite-dimensional manifold — the 
group of volume preserving diffeomorphisms. 

To see this, we have to take the incompressible fluid limit. Of course, a fluid's compress- 
ibility is not a dimensionless quantity. It is a measure of the pressure gradient needed to 
sustain a given density gradient. Since dp/ dp is the (squared) speed of sound, incompressibil- 
ity is not an intrinsic property of the fluid, but rather emerges in the appropriate kinematical 
regime: any fluid behaves as an incompressible one at low enough fluid velocities. In this 
limit, one can consistently restrict to the vortex sector of the theory, i.e. to configurations 
0(x, t) that at fixed time are volume-preserving diffs of x, as we now show. 

5.3.1 Integrating out sound 

It is easy to see that for low speeds and accelerations (i.e. weak time-dependence of 0), 
integrating out classically the compressional modes introduces new Lagrangian terms for the 
vortices that are of order df and higher, whereas the vortex dynamics that we get by simply 
ignoring the compressional modes starts at order d^. Indeed, consider the action for our 
fluid at small velocities. It is more convenient to use the x{(j),t) parameterization for the 
fluid, i.e. to work in comoving coordinates. Moreover, it is instructive to reinstate the speed 
of light c explicitly, to disentangle more transparently relativistic effects from those we are 
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after. We have 

S = -wqC^ j £(pdt det J /((det J~^)^/l - v'^jc^) (63) 

~ Wo y d^4>dt [ - det J /(det J"^) (64) 
+/'(det r^) Wv" + \v^l^\ - \ det J" V"(det r^) v^/c^ ] , (65) 

where J*j is the Jacobian matrix dx"^ /dcjP , and v = dtx{(f),t) is the fluid's velocity field. We 
are giving wq units of a mass density, and / is a dimensionless function of its dimensionless 
argument. The above expression for the action is derived in the Appendix (eq. (86)), or in 
ref. |5]. Now, for nearly incompressible flows, we parameterize the fluid's configuration as 

x{(f),t) = xo{(l),t) + 5x{$,t) , (66) 

where Xq is, at fixed t, a volume-preserving diff of 0, 

{J,y^ ^ ^ , det Jo = 1 , (67) 

and 6x vanishes in the limit of negligible time-dependence of xq. The vorticose motion and 
the fluid's compression are parameterized, respectively, by xo and 6x. At lowest order in 6x, 
this is achieved by imposing that 6x be longitudinal as a function of Xq, that is 

6x = Vo^(fo, t), Vo = ^ = {J^r' ■ . (68) 

OXq 

This corresponds to decomposing the fluid configuration into transverse and longitudinal dis- 
placements, at non-linear order in the transverse ones but at linear order in the longitudinal 
ones. Of course, we may always add background sound waves to our fluid, in which case 6x 
is not necessarily small; but here we are after the compressional modes necessarily associated 
with a given system of vortices by the dynamics. In terms of the above decomposition we 
have 

detJ = l + V^^ + H(V^^)'-(V^V^o^)V--- ' (69) 
and for the velocity field 

v = dt x{$, t) = vo + dtWo'ip ■ (70) 

Notice that these partial time-derivatives are to be evaluated at fixed comoving position 
0, i.e. they are Lagrangian derivatives D/Dt. To integrate out (classically) 5x, we expand 
the action at quadratic order in it. Keeping in mind that will be suppressed by time 
derivatives, the leading terms in a small dt expansion are 



- D 

b ^ + Wq Id Xodt 



lcUvl^?-lcUvl^)^-Vo^-^^vo , (71) 



where 5*0 is the 5a;-independent part of the action, we used /'(I) = 1, /"(I) = c^/c^ (see 
sect. 3|, and we changed coordinates from to xq — for which the Jacobian determinant is 
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one, and upon which we dropped total derivative terms. Notice, crucially, the absence of a 
(Vq^)^ term multiplied by (? rather than by c^, which can be parameterically smaller than 



(? . Naively such a term should be there, from expanding the first term in eq. (65). Yet we 
know that precisely the coefficient of the compressional mode gradient energy defines the 
propagation speed of sound. Therefore terms of order must cancel out, and indeed 



they do — like they did for eq. (10). Moreover, the second term in the action is the only one 
suppressed by c^, and is clearly associated with relativistic effects. For a non-relativistic fluid 
with Cg -C c, it can be safely ignored. On the other hand, for relativistic fluids [cg ~ c) with 
non-relativistic flows [v <C c), such a term should be kept, being of the same order as the 
other ones. Finally, at zeroth order in bx we can replace the Lagrangian derivative D/Dt 
with that associated with the incompressible flow, 

#7 = 1 .+-^o-Vo. (72) 

lyQt Ot const xq 

Varying the action with respect to ip and solving for it, we get 

, 1 1 1^ 1 1 ^ £>o , „„, 
^' = -2Vl? + ^5^^o-:^"«- (73) 

Plugging back into the action (and dropping the subscript '0') we get the effective action for 
the vortex variables, 



S'eff = So + AS (74) 
So = Wo d^(f)dt 



1 2 IcV^ 



AS ~ Wo Sxdt 



2 8c2 8 c 

1 civ* 1 V 1 - D \ 11/1- D ^\ 2 
V ■ —v] + -— — V ■ —V 



c* 2cnV2 Dt ) 2cHV2 Dt 



(75) 
(76) 



where in 5*0 we dropped the rest-mass, w-independent contribution. As advertised, inte- 
grating out the compressional modes introduces 0[df) corrections to the vortex action. At 
small velocities and accelerations we can just restrict our original Lagrangian to the volume- 
preserving configurations, that is we can just use the lowest order term in So- In this 
approximation, and in comoving coordinates, the dynamics are free, but of course we have 
a very non-trivial constraint on the configuration space: x(0, t) cannot leave the space of 
volume preserving diffs. 

Before moving on, it is worth pointing out a couple of features of the computation we 
just performed. The 0[df) correction we computed is local in time but non-local in space. 
Some form of non-locality was to be expected since we are integrating out a massless degree 
of freedom. Yet locality in time follows from our expanding in powers of dt- That is, our 
procedure is equivalent to expanding the sound-wave propagator at low frequencies but finite 
momenta: 

At any finite order such an expansion is analytic in w, and thus local in time. In this sense 
the momentum here is playing the role usually played by the mass when one integrates out 
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a heavy particle in a relativistic QFT. For a given typical frequency w, at momenta of order 
of CsCJ, that is at distances of order l/csW, we expect the expansion to break down. Indeed, 
we know that no matter how slow our vortices, since sound waves are gapless, there will be 
some sound emission. To see (actually to hear) this sound one has to go to the so-called 
wave- zone, which is indeed l/csCJ away from the sources. In a finite volume however, the 
wave-zone may just not be there. Equivalently, there is a gap in the frequency of sound- 
waves. Moreover, even in infinite volume, since interactions involve derivatives, the sound 
emission will be suppressed at very small frequencies and momenta. One can therefore make 
sound emission arbitrarily slow by working at low vortex speeds, thus making the procedure 
of integrating out sound-modes perfectly sensible. It would be interesting to carry out this 
program more systematically, and characterize sound wave-mediated interactions between 
vortices and the sound emission by vorticose motions. We leave this for future work. 



5.3.2 The absence of spontaneous symmetry breaking 

We can now formulate precisely the question of spontaneous symmetry breaking for our fluid. 
We want to compute the amplitude to propagate from x = (p to x = ^{4>), with ^ a volume 
preserving diff: 

(e(0),T|0,O) = / I)fe^^[^] OC e^^[^^'] (78) 

If at late times all volume-preserving diffs get populated with equal probability, there is no 
spontaneous symmetry breaking and our perturbative analysis in the broken phase is not 
applicable. 

Unfortunately, we are not able to carry out this computation for a completely generic 
final configuration ^(0): given the simple free dynamics of the system, we can confidently 
say that there will be a classical solution with time-independent velocity field v{x) evolving 
from X = (f) to X = C,{(f))- However, because of the complexity of the manifold where motion 
takes place, determining the velocity field that connects the initial and final configuration 
would be quite hard. And without knowledge of the velocity field, we cannot compute the 
classical action. 

Nevertheless, to assess the question of spontaneous symmetry breaking we do not need 
to be completely general. We can focus, for instance, on sample vortex configurations with 
a high degree of symmetry. For example, we can consider ^(0) to be a rotation around the 
z-axis, of an angle Acp that depends on the distance r from the axis (dropping to zero above 
some distance R), localized in a region of length L S> -R in the z direction. The classical 
solution connecting the initial configuration to this is a vortex in constant rotation with 
r-dependent angular velocity, 

^ir) = ^ (79) 



6c 



Similarly to the sample f{x) considered in sect. 5.2 the final configuration considered here corresponds 



to exciting an 0(1) number of degrees of freedom. In a region of sizes L, R, and R we are considering 
essentially just the 'fundamental harmonic', with kz ^ l/L and kx ^ ky ^ 1/R. This makes the comment 



at the end of sect. 5.1 immaterial for our purposes. 
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localized in a cylinder of radius R and height L. From eq. (75), the classical action for this 
solution is 

Sci[Aip] c^2n J dzrdrdt ^Wor^co^ir) = j^n J dzdrr^ALp'^{r) ~ Wq ^ , (80) 

where A(p is the typical overall rotation of the final configuration in the region of interest. 
The amplitude we are interested in therefore is 

{^rotation of A(y9 in i?, L ; T |0, O) ~ exp 2 — — — A^p . (81) 

Angles of order one get populated in a cylindrical region of any given size R, and L if we 
wait long enough: 

T > woR'L . (82) 



At late times the symmetry is completely restored. Notice how the estimate (82) matches 
precisely what we could have inferred from the linearized statement (49). 



6 Quantum viscosity? 

Even if the above arguments alleviate the concerns raised by the perturbative analysis of 
sect. [4], we are still left with indications that the effective theory at hand may not be unitary. 

A standard effective field theory can be unitary at low energies thanks to the decoupling 
of the short distance degrees of freedom. That is, thanks to the fact that to excite the 
microscopic degrees of freedom one is neglecting/being agnostic about, one needs a non-zero 
energy. There is a gap, and as long as one works below the gap, the long-distance degrees 
of freedom are sufficient to parameterize the dynamics. 

In our case, we are not in such a good shape. To excite very microscopic vorticose 
deformations of the fluid we need no energy at all. High momenta are not associated with 



high energies. This is also evident from eq. (49 ), where energy and momentum play strikingly 



different roles. This suggests that for any given cutoff in momentum space, the effective 
theory may not be exactly unitary. No matter how small the energy of the process under 
consideration, there may be a non-trivial probability flow across the cutoff. From the effective 
theory viewpoint this should look like dissipation. In fact, for classical turbulence in a viscous 
fluid this is exactly what happens — viscosity drives vorticity from large scales to smaller and 
smaller ones, down to the UV cutoff of the fluid description (the mean free path of the 
underlying microscopic system.) Could it be that in our case we have some sort of quantum 
contribution to viscosity due to this non-decoupling of micro-vortices? How can we test this 
conjecture? 



7 Discussion and Outlook 

Our flndings raise more questions than they answer. The perturbative analysis about the 
naive, semiclassical vacuum 0^ = indicates that the ordinary fluid effective fleld theory 
is strongly coupled at all scales, and thus inconsistent (sect. |4]). On the other hand, a more 
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careful non-perturbative study of the theory's quantum-mechanical vacuum shows that this 
has essentially nothing to do with the semiclassical one, and suggests that the naive per- 
turbative degrees of freedom and their dynamics have no quantum-mechanical counterpart 
(sect. [s]). One is thus tempted to welcome the latter conclusion and drop the former, and 
simply ignore the perturbative results. In particular, the strong coupling problem we isolated 
in sect. |4] might not be there — being associated with excitations that might themselves not 
be there. Also, in the x{$, t) parameterization of sect, [s , the vortex dynamics are essentially 
free — the only thing resembling an interaction is the volume-preserving constraint. It is not 
clear what 'strong coupling' would mean in such a description of the system. 

However, there is a number of confusing aspects that suggest that this optimistic attitude 
may be naive. The first is that even for theories where Coleman's theorem applies, one 
can recover correct information about physical quantities — like the spectrum for instance — 
by doing perturbative computations about the wrong semiclassical vacuum [9] (see also a 
related discussion in \IQ\). We do not know yet whether, and to what extent, the results of 
[9] apply in our case. If they do, they would demote our Coleman theorem-like result to a 
somewhat formal statement about the vacuum structure in the far infrared, with less crucial 
consequences for more local physics. 

The second confusing fact — which points in the same general direction as the previous 
one — is that for tiny but finite ct we have a well defined effective theory that is perturbative 
up to some finite energy/momentum scale. For this effective theory there is no funny Coleman 
theorem-like behavior. All the symmetries that are spontaneously broken classically, remain 
so quantum-mechanically. The associated Goldstone bosons are described precisely by our 
perturbative analysis of sect. |4j On the other hand, as we stressed above, if ct <^ cs we 
expect this theory to be physically equivalent to the ct = one for local questions, which 
do not rely crucially on the precise nature of the asymptotic states of the theory. While this 
kind of logic may be misleading for massive gauge theories or massive gravity, where having 
or not having the mass really determines the number of local (as opposed to asymptotic) 
physical degrees of freedom, here there is no such subtlety. The transverse degrees of freedom 
are perfectly physical even for vanishing ct', they have non-vanishing conjugate momenta and 
are thus standard Hamiltonian degrees of freedom. Only, they do not feature wave solutions. 

It thus seems that for local questions the perturbative analysis of sect. |4] should be 
perfectly fine. Of course the S'-matrix is not a local quantity, and it may well be that the 
pathologies we encountered there are irrelevant for local questions. But we find it unlikely: 
if an effective theory exhibits strong coupling in scattering processes at some energy and 
momentum scales, it is probably useless for computing local correlation functions at the 
corresponding length scales. To settle the question one should compute directly local n- 
point functions, and see whether the perturbative expansion breaks down there. However, 
for the S'-matrix we have a very powerful property — unitarity — that makes the tree level 
sufficient for such a question. For local correlation functions instead, one should really 
ascertain the validity of the perturbative expansion by computing loop corrections — which, 
given that Lorentz invariance is spontaneously broken, is certainly doable but somewhat less 
transparent than usual. We leave this to future work. 

Our results invite one to focus on correlators of quantities that are invariant under volume- 
preserving diffs. For instance, for a free massless scalar in 1 -|- 1 dimensions, even though 
Coleman theorem applies, correlators of shift-invariant quantities are perfectly well-defined. 



28 



and match what one would naively expect by doing perturbation theory about the wrong 
classical vacuum where the shift symmetry is spontaneously broken. Whether this property 
survives the inclusion of interactions and the generalization to our fluid, we do not know. 
If all diff-invariant correlators are well defined for our fluid, then one could in principle 
decide that those are the only observable quantities. This would correspond to gauging 
the problematic volume-preserving symmetry. This is certainly an interesting possibility to 
consider, and maybe the theory defined this way would be consistent. Yet it is not clear 
to us what resemblance it would bear with a physical fluid: as we tried to make clear, at 
the classical level the volume preserving diffs are not a gauge redundancy — they are real 
symmetries acting on physical and measurable degrees of freedom. 

On a different note, even if the effective theory is not strongly coupled, we find it inter- 
esting that microscopic vortices may impair its unitarity (sect. |6|. We plan to make this 
statement more systematic and quantitative. In particular, it would be interesting to un- 
derstand to what extent this effect can be parameterized as a new contribution to viscosity, 
and whether it has any relation to the conjectured viscosity-over-entropy bound [TT] . 
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A Expanding the Lagrangian 



To carry out the expansion of the Lagrangian in fluctuations, eq. (16) is a slightly more 
convenient starting point than eq. ([T]). More importantly, it is particularly convenient to 
'pull out' the time derivatives from B: 

B = detd^(l)jd''(l)j 
det ^^-^^ 



= det (d(f)'^ ■ (1 - {id(/)^r' ■ i) ® ((90^)"' • ■ d(j?^ 
= (det (90)2 det (1 , (83) 

where we defined the matrix {d(f))ij = di(j)j, and the vector 

V = {d(f)'^)-^ ■ $ (84) 
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Notice that v is the usual fluid velocity field — hence the name. The last term in (83) is easy 
to compute, for instance by going to a basis where, locally, the x axis is aligned with v. We 
get 

B = {detd(j)y{l-\v\'^) (85) 

and therefore 

£ = -wo/(det90 - (86) 

Notice that we are essentially reproducing eq. (88) of [5j. We now have to expand / in 
powers of its argument, and its argument in powers of vr. The benefit of doing the expansion 
this way is that one now only has to expand in vr functions of (90 = 1 + dn, rather than of 
dn(j)d'^(j) = 1 + dn + dn^ + dn^dn + tttt. This simplifies the algebra considerably. We just 
need 

{dcj))-^ - l-dTT + dTT^ (87) 
det50 = 1 + [dn] + ^{[dn]^ - [dn^]) + l{[d7rf - 3[dTT][d7i''] + 2[d7r^]) (88) 

(the determinant of an n x n matrix stops at n-th order.) A straightforward Taylor-expansion 



of (86) up to fourth order then yields eq. (17). 



It is worth pointing out two sources of non-trivial cancellations, with important physical 
consequences. The first forces all Lagrangian terms that do not involve time-derivatives to be 



weighed by = /"(I) or by higher derivatives of /, as manifest in eq. (17). That is, no such 



term is coming from the expansion of (86) at first order in /'s argument, i.e. with a coefficient 
/'(I) = 1. The reason is simple: neglecting time-derivatives, the term proportional to /'(I) 
would be 

C D -Wof'{l) ( det 90 - 1) , (89) 

which is a total derivative: 

detd(j) = eed(f)d(f)d(f) = d{ee(j)d(j)d(j)) (90) 

Since one expects higher derivatives of / to be naturally of order c^, for a non-relativistic fiuid 
this cancellation has the effect of weakening the interactions considerably, or equivalently 
of raising the strong-coupling scale compared to what one may have naively guessed before 
carrying out the expansion. 

The second cancellation involves the transverse phonons only, and has also the effect 
of weakening some interactions and correspodingly raising the vortex strong-coupling scale. 
Consider an interaction term with spatial derivatives only, and assume that at least one 
of the phonons entering the corresponding vertex is transverse. As not manifest from the 



Lagrangian (17), such a vertex yields zero. The reason is that we can perform a non-linear 
field redefinition that makes vortices disappear from all Lagrangian terms without time- 
derivatives. The trick is to define vf so that 

det 90 = 1 + V ■ 7f (91) 



exactly. That this is possible follows from eq. (90) — we may as well call the total derivative 



on the r.h.s. 1 -|- V • vf. This matches our original definition of vr at linear order, and as a 
consequence it does not affect the S'-matrix. But now it is clear that vortex interactions will 
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only come from the part of (86), and will thus involve at least two time-derivatives. The 
downside is that in these variables the structure of the Lagrangian will be more complicated 
than eq. (17); in particular we will not have exactly one derivative acting on each field. For 
this reason we stick to the original definition of the phonon field and to eq. (17), but we should 
expect non-trivial cancellations when computing S'-matrix elements involving vortices, as we 
indeed find in sect. HI 



B The /S-matrix, cross-sections, and decay rates 

Here we briefly review the standard relativistic formulae for the S-matrix and related phys- 
ical quantities like cross sections and decay rates, and derive the modifications needed for 
applying them to our c 7^ 1 case. The rules we will derive are straightforwardly generalizable 
to the case of different fields with different propagation speeds. 

We will borrow the conventions of Peskin-Shroeder [7] . In particular, we use the so-called 
relativistic normalization for one-particle states: 

{p\q) = {2E) {2^)H\p- q) , (92) 

(we are suppressing spin labels — their inclusion is straightforward) and of course the vacuum 
state |0) is normalized to one. This way, a relativistic canonically normalized scalar field 
0(x) obeys 

(0|</.(x)|p) = e-^(^*-P^-^^), (93) 

and consequently the momentum-space Feynman rules assign a wavefunction one to external 
spin-0 states. One thus has that for a 2 — )■ rij scattering process the infinitesimal cross section 
is 

= rr^ri^u. V~\ \^AB^f\^ dli^^. . (94) 

Here M.i^f is the amplitude computed according to the standard relativistic Feynman rules, 
and defined by 

(gi . . . \ {S - l)\pi . ..pn,) = (27r)^ 5'^ (momentum) ^(energy) • iMi^f , (95) 
and the dUnj^ is the relativistic final-state phase-space: 

rfH„^ = (27r)^ 5^ (momentum) 5 {energy) ■ ^ JJ 2^) ^^^^ 

Finally, |f a — "^^bI is the relative velocity between the two colliding beams as measured in the 
lab frame. Likewise for a 1 — i- decay process, the infinitesimal rate is 

dr = -]-\MA^f\'dUr,^. (97) 

First, let us check the dimensions of these quantities, by keeping h dimensionless but the 
speed of light dimensionful. That is, let's give energy = 1/time and momentum = 1/length 
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different units. From their definitions, eqs. (95 96), for the amphtude and phase-space 
element we get 



[M] = Ek\E/k^)^' , [dU] = {Ek'^)-\E/k'^)"'' . (98) 

The cross-section and decay rate thus have dimensions 

[da] = k-^ = area , [dV] = E = 1/time , (99) 

as they should. This means that the above fomulae are aheady dimensionally correct with 
no need of exphcit powers of the speed of hght. 

Next, we notice that nowhere is Lorentz invar iance assumed in deriving the Feynman 
rules and the above expressions for a and F. This is evident e.g. in the derivation of ref. [7], 
apart from the relative velocity factor in the cross section. However that too is independent 
of Lorentz invariance, for it arises from the integral of an energy delta-function over the 
longitudinal (w.r.t. to the collision direction) momenta of the incoming wave-packets: 

j dk\dkl6(k\ + k'B-P})S(EA + EB-E,) = j dk\S(EA + EB-E,)\^,^^p._ 

-\ ' ' (101) 



For each wave-packet, the derivative of the energy w.r.t. the corresponding momentum is the 
wave-packet's group- velocity, independently of the actual form of the dispersion law E{k). 
The above thus yields the factor l/|fA — vb\ in the cross section, regardless of Lorentz- 
invariance. 

The bottom line is, much ado about nothing. We can use the standard relativistic 
Feynman rules and formulae for infinitesimal cross-sections and rates for our non-relativistic 
case as well, with no modifications, even when different fields have different speeds. The only 



subtlety we should keep in mind is that canonically normalized fields obey eq. (93), times 
possible polarization factors for non-scalar particles. This means that a scalar field thus 
normalized should appear in the action as 



S = j d'xdt + . . . , (102) 
so that single-particle states are eigenstates of the free Hamiltonian with the right energy: 

Mp) = [J d'xl4>' + ---) \P) = E{p)\p) . (103) 

As a check that these conclusions make sense, we estimate the cross section for sound 
wave-sound wave elastic scattering and show that, indeed, we have strong-coupling at the 



correct energy. From Feynman rules applied to the Lagrangian (17) we have 

Mr^cl— (104) 

Wo 
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where the factor of Wq comes from the non-canonical normahzation of vr^. The final state 



phase space (96) is of order 

^f-V^^i^'/^T^ (105) 



and the relative velocity is of course 2cs, so that the cross-section (94) is 



) • (106) 

fc^ \WoCsJ 

This agrees with our estimates of sect. |3.1) — see the last paragraph of sect. 4A_ 



C Phase space 

We are mostly interested in a two-particle final state, possibly with two independent propa- 
gation speeds. The infinitesimal phase space is 

dn, = (2.r 5HP - ,1 - S^E-E,-E,).^,^,±-±-, (107) 

where E and P are the total energy and momentum. The integral in q2 eliminates the 
momentum-conservation delta-function. Then we are left with 

''n.^^,;^;,,^*!^ (108) 

with the understanding that E2 be evaluated at q2 = P — qi- We have 

6{E -E,- E2) = (109) 

I dqi (9(32 dqi \ 



and 

dq2 _ d\P - qi\ qi-Pcos9 
dqi dqi q2 



:iio) 



where 6 is the angle between qi and P. On the other hand, the derivatives of the energies 
w.r.t. the corresponding momenta are the particles' group velocities. Integrating over qi we 
thus get 

dll2 = —,4^, r (111) 

IGn^ E1E2 \c1q2 + C2qi — C2Pcos^^| 
For a linear dispersion law like in our case, Ea = Caqa, we finally have 



dQ 1 qi 
IGn"^ C1C2 \c1q2 + C2qi — C2Pcos6'| 



d^2 = — -2—7—^^-^7— 7^- (112) 



In special circumstances there are further simplifications: 
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i) For scattering processes at zero total momentum, we can set P = and qi = q2- We 
get 

dU, = ^ r^——^ iP = 0). (113) 

167r^ ciC2(ci + C2) 

ii) For decay processes at finite total P, but when one of the final particles is much slower 
that the other, barring an hierarchy between qi and q2 we have 



dn^ ^ ^ ^ (C2 « ci) . (114) 

IbTT^ C1C2 q2 



Of course the ratio qi/q2 depends non-trivially on the angle 9 we are supposed to 
integrate over — which we can take to be the angle between g! and P. We have: 

1-2^ fe«c.P^O). (115) 

Overall we thus get 

'^n^-^^S^ ta«c,P/0). (116) 
Notice that this is regular at 6' = 0, thus making our 'barring an hierarchy . . . ' approx- 



imation under control. That is, eq. (116) is the correct phase-space element at lowest 
order in C2/C1. 
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